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Correction  for  Ionospheric  Refraction  for  COBRA  DANE 


1.  INTRODUCTION 

Refraction  effects  within  the  ionosphere  degrade  measurements  of  range  and 
bearing.  For  the  COBRA  DANE  radar  these  effects  create  errors  that  exceed  the 
required  metric  accuracy;  t!iSr  .‘fore,  a correction  must  be  made  for  each  target. 

In  developing  the  correction  procedure,  consideration  had  to  be  given  to  the  present 
and  future  requirements  of  the  user,  the  constraints  within  the  particular  radar, 
the  forecast  capabilities  of  the  Air  Force  Global  Weather  Central  (AFGWC),  the 
Inherent  limitations  of  the  ionosphere,  and  foreseeable  techniques  that  will  des- 
cribe the  ionosphere  present  at  the  radar.  A simple  method  of  correction,  based 
on  three -element  vectors,  provides  the  equivalent  of  a look-up  table  for  the  entire 
surveillance  volume  of  the  radar.  This  correction  procedure  satisfies  the  present 
system  requirements  and  has  flexibility  for  expansion  as  techniques  improve  or 
requirements  change. 

Once  a month  the  vector  model  is  calculated  off  line  by  the  AFGWC.  Used 
alone  by  the  COBRA  DANE  processor,  it  removes  about  75  percent  of  the  monthly 
mean  ionospheric  effect.  Since  the  requirements  for  metric  accuracy  on  single 
missions  will  become  more  critical  for  the  enhanced  ionosphere  at  the  next  solar 
maximum,  provision  has  been  made  for  refining  the  correction.  This  is  now  done 
in  two  ways;  first,  from  worldwide  ionospheric  observations,  the  AFGWC  deter- 
mines and  supplies  a dally  update  factor  to  scale  the  correction  hour -by-hour  so 
(Received  for  publication  18  November  1077) 
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as  to  more  accurately  represent  the  current  ionosphere;  second,  using  a dual  fre- 
quency range  measurement  to  certain  satellites,  the  radar  determines  a scale 
factor  that  provides  a more  precise  correction  for  a shorter  time  span.  The  accu- 
racy of  this  procedure  with  such  scaling  factors  will  depend  on  the  frequency  of 
update.  If  a scaling  factor  can  be  determined  in  the  area  of  interest  as  frequently 
as  every  15  min,  then  perhaps  95  percent  of  the  ionospheric  effect  can  be  removed. 


2.  GENERAL  CONSIDERATIONS 

The  final  correction  procedure  that  was  delivered  and  installed  in  the  COBRA 
DANE  radar  as  software  for  a real-time  hit -to-hit  correction  for  ionospheric  re- 
fractive effects  is  an  adaptation  of  the  model  designed  to  provide  the  basic  monthly 
mean  correction  for  the  SPADATS  mission.  ^ 

The  rationale  for  the  procedure  has  not  been  changed.  The  customer  require- 
ments for  metric  accuracy  during  the  1980's  could  not  be  satisfied  by  correcting 
for  just  the  expected  monthly  median  refraction  effects,  which  would  have  removed 
only  about  7 5 percent  of  the  day-to-day  variability.  A procedure  was  needed  to 
improve  upon  this  basic  estimate  by  some  update  that  could  follow  local  features  of 
the  current  ionosphere.  Yet  the  constraints  within  the  radar  processor  on  core 
storage  and  computation  time  would  not  permit  consideration  of  any  procedure  that 
was  more  complicated  than  a look-up  table.  Any  mapping  techniques  or  line  of 
sight  ray  tracing  methods  required  real-time  evaluation  of  analytic  functions  that 
so  reduced  throughput  within  the  processor  that  a correction  could  not  be  made  on 
a real-time  hit -to-hit  basis. 

As  a specific  example,  a simple  mapping  algorithm  was  developed  to  represent 
the  expected  median  correction  for  the  ionospheric  range  component  over  the  whole 
field  of  view.  Using  11  Kourler  sine  terms  to  represent  the  longitudinal  variation 
and  11  Tchebyshev  polynomials  to  represent  the  latitudinal  variations,  a correction 
could  be  calculated  that  would  match  the  expected  value  calculated  from  the  original 
data  field  with  an  accuracy  near  5 percent.  Yet  when  this  minimum  slsied  algorltlun 
was  tested  on  the  equivalent  ClXl  radar  computer,  the  computation  took  so  long  that 
a hit -to -hit  correction  could  not  be  made. 

When  the  only  solution  to  a hit -to-hit  correction  la  a quick  access  look-up  table, 
then  a now  problon;  Is  created.  If  a procedure  is  to  bo  developed  that  will  meet 
present  goals  and  grow  as  needs  and  techniques  develop,  the  cell  size  in  the  radar 
volume  must  be  fairly  small.  To  attain  the  goal  of  removing  95  percent  of  the 
variability,  local  features  in  the  ionosphere  have  to  be  introduced  by  some  updating 

1.  Allen,  R.S.  (1974)  An  Ionospheric  Correction  Model  for  the  COBRA  DANE  Radar, 
AFCRL  Memo  i,m77f-2543. 
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technique.  The  look-up  table  procedure  then  has  to  be  quite  fine  in  azimuth,  eleva- 
tion angle,  and  height  of  target— so  fine,  in  fact,  that  it  becomes  unreasonable  in 
terms  of  storage  area  within  the  computer. 

The  proposed  solution  is  an  engineering  com.promise.  The  goal  is  the  metric 
accuracy  that  the  user  wants  for  single  missions  in  the  ionosphere  of  the  early 
1980's.  The  correction  procedure  installed  in  COBRA  DANE  is  a three-vector 
model  that  provides  the  equivalent  to  a fine-scale  look-up  table.  This  vector  model 
is  computed  for  each  month  from  forecasts  made  by  the  AFGWC.  It  is  presently 
in  use  at  the  COBRA  DANE  radar  and  should  be  updated  by  a scaling  factor  to 
represent  overall  ionospheric  conditions  on  the  day  of  operation.  As  techniques 
are  developed  for  scaling  in  space  and  time  to  represent  the  actual  local  state  of 
the  ionosphere  during  missions,  they  should  be  adopted  by  the  radar.  The  basic 
monthly  mean  forecast  is  the  default  correction  when  any  of  the  updates  are  not 
available.  Most  important,  the  software  installed  at  the  radar  is  small  in  total 
storage  words  and  fast  in  computation  time  so  that  a hit -to-hit  correction  can  be 
made. 


3.  INTUITIVE  DEVELOPMENT  OF  THE  VECTOR  MODEL 


The  problem  of  devising  a method  of  approximating  a very  fine  scale  look-up 
table  is  divided  into  two  parts:  the  estimation  of  the  geometric  effect  of  viewing 
targets  along  slant  paths,  and  the  estimation  of  the  ionosphere.  To  meet  accuracy 
requirements,  the  beat  correction  is  made  in  th.at  portion  of  the  field  of  view  where 
it  is  most  needed  and  some  degradation  is  allowed  in  other  .areas.  Since  both  the 
intelligence  and  the  SPADATS  missioms  at  the  COBRA  DANE  radar  require  the 
highest  accuracy  for  targets  embedded  in  the  ionosphere  at  low  elevation  angles, 
the  same  tnodel  can  bo  used  for  both  missions. 


The  ionospheric  and  geometric  effects  cannot  be  considered  separately,  simply 
because  the  ionosphere  has  an  arbitrary  distribution  along  the  slant  path  to  the 
target.  However,  the  day-to-day  variability  about  the  forecast  of  the  monthly  moan 
ionosphere  is  about  10  to  20  percent  for  each  parameter,  much  larger  than  would 
be  expected  from  the  variation  along  the  slant  path.  With  this  justification,  the 
ionosphere  Is  defined  at  the  central  position  of  the  slant  path  and  geometric  effects 
are  calculated  as  if  the  ionosphere  wore  spherically  Gratified. 

The  elements  of  a look-up  table  to  provide  a correction  over  the  entire  radar 
volume  can  be  ordered  by  oidhogonal  unit  vectors  for  azimuth,  elevation  angle,  and 
target  height.  In  this  model  the  ionosphere  Is  determined  near  the  horizon  where 
the  greatest  accuracy  is  desired,  and  the  geometric  effects  are  determined  along 
these  vector  directions  so  that  the  refraction  correction,  C. ,,,  can  be  computed 

ijK 

within  the  radar  by  simple  multiplication: 
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where 

(1)  is  the  magnitude  of  the  correction  in  feet,  at  the 
prime  reference  point,  azimuth  319°  T.  height  1000  km, 
and  elevation  5.3°. 

(2)  The  obliquity  factor,  S^,  is  a 33-element  vector  for 
the  variation  with  elevation  angle  evaluated  at  the 
reference  azim'jth  and  height  (310*  T and  1000  km), 
normalized  to  1.0  at  5.3°  elevation, 

(3)  The  partition  factor,  Y^,  is  a 33-element  vector  for 
the  variation  with  the  height  of  the  target  evaluated 
at  reference  azimuth  and  elevation  angle  (319°  T and 
5.  3°),  normalized  to  1.  0 at  1000  kni. 

(4)  The  ionospheric  factor.  is  a 33 -element  vector  for 
the  variation  with  azimuth  evaluated  for  the  vertical 
ionosphere  at  the  center  of  each  azimuthal  sector  near 
the  horizon  (5.  3*  elevation,  420  km  height).  This  is 
normHU..cd  to  1.0  at  boicsight,  SIO”  T,  and  has  been 
weighted  for  the  variation  of  S (5.  3*)  wltl\  azimuth. 

(5)  The  update  factor,  U,  represents  ll)e  cnirrent  ioiwsphere, 
either  the  scale  factor  from  the  AEGWC  or  a recent  local 
measurement. 

From  the  101  stored  values  (C  + S i Y ♦ E * U'.  the  radar  computer  can 

o 

generate  the  equivalent  clomontB  of  a 35,937  value  look-up  table.  As  a result  of 
the  definition  of  tlio  reference  conditions,  the  correction  is  most  accurate  in  the 
precision  sector,  and  for  those  targets  that  require  the  greatest  correction.  As 
the  target  coordinates  shift  away  from  the  reference  point,  the  vectors  depart 
froni  their  true  values,  but  the  magnitude  of  the  correction  decreases  so  that  a 
larger  percentage  error  can  still  be  tolerated.  For  instance,  for  a target  near 
80*  elevation,  the  ionosphere  may  be  2 hou>‘a  renmved  from  the  reference,  but  the 
obliquity  olcmerst,  S,  at  80*  is  about  one-third  that  at  5.  3*. 

The  result  is  a very  small,  very  rapid  computer  program  for  generating  hit- 
to-hit  corrections  with  fine-scale  co\’eragc  over  the  entire  radar  volume  and  an 
error  that  is  minimized  at  the  maximum  correction.  The  value  in  any  cell  of  the 
equivalent  look-up  table  is  expected  to  be  within  a few  percent  of  ti>at  expected 
when  U8lt>g  the  monthly  mean  prediction.  In  addition,  the  correction  for  refractive 
effects  can  be  updated  for  local  conditions. 


8 


4.  COMPUTATION  OF  VECTOR  ELEMENTS  BY  APGWC 


A Fortran  subroutine  (subroutine  DANE)  has  been  developed  which  computes 
the  required  vector  elements  to  be  used  at  the  radar.  The  subroutine  is  tailored 
to  COBRA  DANE  by  internal  Fortran  data  cards  that  specify  the  33  ionospheric 
reference  points  at  the  center  of  each  a2imuth  sector  and  the  geometric  values 
pertinent  to  the  33  elements  of  the  obliquity  and  partition  factors.  It  is  initialized 
at  the  AFGWC  by  specifying  the  desired  month  and  the  estimated  12 -month  running 
mean  of  the  sunsfiot  number  and  the  2800  MHz  solar  flux.  The  subroutine  then  calls 
a standard  subroutine  available  at  AFGWC.  such  as  the  Damon /Hartrampf  or  Bent 
procedures,  for  building  the  expected  electron  density  profile  at  the  reference 
{joints.  This  standard  subroutine  calculates  and  retai'os  the  vertical  total  electron 
content  (TEC  in  el/cc)  up  to  1000  km  and  the  height  (hmF2  in  km)  of  the  peak  of  the 
F2  region. 

Subroutine  DANE  then  uses  hmF2  at  boreslght.  sector  17,  to  calculate  the 
obliquity  factor,  SHAH: 

SHAH  (EL)  « (2.02-0.  0053(hmF2-290)| 


for  each  of  the  33  stored  values  of  elevation  angle  (EL).  .ALso,  for  boreslght,  the 
partition  factor,  Y,  is  calculated  for  each  of  the  stored  values  of  height  (HT): 


Y(HT) 


1. 


-B/ 


H + 


/HT-100. 
llmiPI-lDO  I 


where 


B • 1.81  T 0.001  {hmF2-290)  . 


Note  that,  by  definition,  Y approaches  1.0  at  1000  km.  The  reference  value,  C^, 
is  then  calculated  for  each  hour: 


C » 
o 


1.3225X10 

p 

op 


-10 


TEC{n)SHAH{5.3*)  (in  feet)  . 


whore  the  raciar  operating  frequency,  nominally  1275  MHa.  The  obliquity 

factor  is  then  nomialUcd  to  l.O  at  the  reference  elevation.  5.  3*. 


S(EL)  = SHAH(E)/SHAH(5.3*)  . 


and  the  TEC  at  each  azimuthal  position,  (AZ).  Is  normalized  to  boreslght  and 
weighted  by  the  obliquity  factor  for  each  azimuth  at  5.  3*  elevation. 


_ TEC(AZ)SHAH(AZ.  EL) 

° TEC|17)SHAH(A'Z.  5:3^7  * 


to  provide  the  33  elements  of  E,  the  azimuthal  vector,  which  is  normalized  to  1.  0 
at  boresight. 

The  output  of  subroutine  DANE  is  a set  of  10  punched  cards  for  each  integer 
hour  of  Universal  Time  (0-23),  as  shown  in  Figure  1. 
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Figure  1.  Format  of  Cards  PrtHiuced  by  GWC/AFtlL  Frogram 


The  normalising  condition  for  the  elevation  angle  and  for  the  ionospheric 
control  point  wan  clK>son  st  5.  3'  to  provide  an  accurate  correction  in  the  region  of 
greatest  need  and  to  diffuse  the  ev-or  in  the  0 to  15  deg  ele\'at(on  region.  The 
analytic  forma  of  she  olilicjulty  and  partition  factors  vrere  derived  from  an  AFCU, 
study  of  electron  density  profiles  obtained  fi'tJm  the  Millstone  Hill  incoherent  scat- 
ter radar.  These  functions  ore  a best -fit  to  the  exjjected  <loytime  values  of  a high 
mid -latitude  ionosphere. 


5.  rxtMRn  AnoN  of  refraction  correction  by  cobra  dank 

A pre'ets  of  the  subroutine  used  by  the  COBRA  DANE  radar  to  correct  for  the 

iortosphcrl'  component  in  the  range  measurement  le  as  follow#: 

(1)  Each  hour  the  central  processor  stores  101  word# 

(C  ♦ S * Y'  ♦ E * U)  to  provide  vector#  for  the  refraction 
0 

correction,  supplied  on  card#  for  each  month  by  AFGWC. 
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(2)  The  hourly  value  of  the  update  factor  supplied  via  Autodin 
by  the  iiFGWC  is  resident  in  storage.  When  a pulse-pair 
measurement  has  been  made  with  one  of  the  calibration 
satellites,  it  is  automatically  entered  into  the  subroutine 
and  used  in  place  of  the  AFGWC  value.  When  the  time 

since  e calibration  exceeds  3 hours,  or  if  the  user  determines 
it  is  erroneous,  the  update  reverts  back  to  the  AFGWC  value. 

(3)  For  each  hit  or  batch  of  hits,  the  subroutine  is  entered  with 
the  readily  accessible  coordinates  of  target  position:  the  sine 
of  the  elevation  angle,  the  range,  and  the  tangent  of  the  angle 
off  boresight.  For  the  refraction  correction  these  are  re- 
defined as  follows: 

• An  index  for  the  elements  of  the  vector  S is  defined 
by  dividing  sine  E to  provide  33  equal  increments  of 
sine  E from  0*  to  80*: 

S - slnE/0.Q3077  + 1. 

In  addition  to  being  a fast  computation,  this  also 

provides  greater  definition  near  the  horlton  whore 

It  is  most  needed. 

^ •) 

9 As  prot>o8ed  by  Kate  and  Boak.  “ it  is  convenient  to 

order  the  variation  with  height  by  the  transformation 

I • K{(r/H)^  + 3(r/R)sii»eE)  . 

The  constant  K Is  defined  when  I is  33  at  1000  km  and  H 
is  the  radius  of  the  earth.  This  transformation,  nearly 
linear  in  helgrrt  intet-val.  avoids  transcendental  functions. 

The  average  .spacing,  ai.out  31  km.  is  about  half  of  the 
char,ticterl8tic  scale  slae  near  the  twsk  of  the  F region, 
so  again  the  . ell  siac  is  fine  scale. 

• Aisother  easily  computed  index  pi^oposel  by  Kata  and  Hoak*  is 

e ' tanAAZ/d.  *^0.2  lan^AAZ)  . 

where  tanAAZ  is  the  tangent  of  the  angle  off  boresight. 

This  provides  a fine  spacing  of  al>oul  4*  of  aairnuth  in  (he  precision 

2,  Katz,  A.H.  , aivd  Boak.  Ill,  T,  I.S.  (1P771  Interface  Control  I>ocu!T!cnt  for  the 
COBRA  DANE  Ionospheric  Hansje  Con-rclion 
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sector,  increasing  to  about  5“  at  the  edge  of  the  coverage 
area. 

A The  correct  elements  of  the  vectors  S,  Y,  and  E are 
calculated  by  linear  interpolation  between  the  stored 
values.  The  range  correction  is  calculated  by  simple 
multiplication: 

C = C^S(3)Y(I)E(e)  . 
and  is  returned  by  the  subroutine. 


6.  UPDATE  FOR  CURRENT  IONOSPHERIC  CONDHTONS 

At  present  there  are  two  update  procedures,  both  of  which  are  a scaling  of  all 
of  the  corrections  in  the  field  of  view  by  the  single  update  factor  U.  The  basic 
default  update  is  that  supplied  by  APGWC  from  their  world  observations.  This 
update  factor  is  a set  of  hourly  factors,  provided  on  a daily  basis.  It  is  defined  as 

J (HR) = ETEC/PTEC  , 

where,  for  each  hour  of  the  predicted  period,  ETEC  is  the  new  estimation  and 
PTEC  is  the  baseline  prediction  supplied  by  the  card  to  the  radar,  both  for  the  ver- 
tical electron  content  at  the  reference  coordinates  of  the  precision  sector,  63.  3°N. , 
150,  The  24  values  of  the  update  factors  are  supplied  to  the  COBRA  DANE 

radar  via  Autodin  in  decima'  notation  and  are  entered  directly  into  the  computer 
by  the  oper.ator. 

A short-time- span  update  is  derived  by  the  radar  from  pulse-pair  measure- 
ments of  the  range  to  certain  satellites.  This  can  be  made  any  time  the  radar  can 
see  any  one  of  four  calibration  satellites,  each  with  several  passes  per  day,  in  the 
sector  299°  to  337  °T  in  azimuth  and  0°  to  42°  in  elevation.  A procedure  for  doing 

this  has  been  worked  out,  tested,  and  Installed  in  the  software  at  the  radar  by  the 
2 3 

contractor.  In-house  studies  at  AFGL  have  shown  that  this  update  provides  a 
very  significant  reduction  in  the  range  residuals  over  short  time  intervals,  but 
that  Ihe  factor  should  revert  back  to  the  current  AFGWC  update  when  the  pulse- 
pair  update  is  older  than  3 hours.  In  practice,  the  radar  tries  to  renew  the  pulse- 
pair  update  within  the  limits  of  1/2  to  3 hours. 


3.  Donateili,  D.  E.  (1977)  Reduction  of  the  Uncertainty  of  Radar  Range  Correction, 
AFGL -TR -77 -0125. 
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7.  EFFECTIVENESS  OF  UPDATING  TECHNIQUES 


A study  of  the  temporal  variability  of  the  total  electron  content  (TEC)  of  the 

3 

ionosphere,  as  measured  from  a mid-latitude  station,  has  yielded  results  that  can  be 

applied  directly  to  the  COBRA  DANE  radar  to  estimate  the  effectiveness  of  the 

updating  technique  using  pulse -pair  measurements.  Range  correction  is  directly 

proport<orc.i  to  TEC,  and  for  an  L-band  radar  the  proportionality  constant  for 

1 R 9 

range  correctioi  , infeet,  is  0.7  X TEC  in  units  of  10  el/m  . 

The  average  day-to  -day  variability  of  TEC  is  20  to  25  percent  of  its  monthly 
median  in  d£,y'ime  and  30  to  33  percent  at  night.  This  percentage  error  is  nearly 
Independent  of  seasonal  and  solar- .ycle  variations.  The  absolute  error  is  2 to  3 
times  greater  in  daytime  than  at  nighttime. 

An  updating  technique  combining  the  use  of  pulse-pair  measurements  with  a 
model  prediction  of  median  range  correction  can  reduce  the  residual  error  sig- 
nifh  tt-  - particularly  in  daytime  when  the  error  is  greatest.  Evaluation  of  an 
upd.-.inj  technique  using  9 years  of  data  from  Hamilton,  Massachusetts  (in  a 
simulation  of  actual  use  by  a 1300  MHz  radar  on  a target  at  1000  km  altitude,  5 deg 
elevation  angle),  has  shown  that  an  updated  prediction  of  the  median  range  correc- 
tion can  reduce  the  residual  error  by  60  percent,  even  after  1 hour.  In  the  daytime, 
when  the  largest  correction  Is  needed,  the  residual  error  can  he  reduced  30  percent 
or  better,  even  after  3 hours.  Attempts  to  use  the  updating  factor  for  a period 
longer  than  3 hours  results  In  an  error  that  is  of  the  same  order  as  the  error  that 
results  by  simply  using  the  predicted  tnedlan.  Near  sunrise  and  sunset  and  during 
severe  magnetic  disturbances,  all  of  which  are  periods  of  rapid  change  in  TEC,  the 
same  degree  of  error  reduction  can  be  maintained  by  reducing  the  Interval  for  deter- 
mining new  updating  factors  to  about  15  to  30  minutes. 

The  results  of  this  study  are  summarized  in  Table  1 for  the  solar  maximum 

condition.  155,  K allO.  and  the  solar  nunimum  conditions,  "S*?  1,  R ’10,  where 
z z 

S is  the  12-mcnth  running  mean  solar  flux  at  2800  MHs,  and  H la  the  12-month 

' z 

running  mean  sunspot  nvimber.  The  values,  In  feet,  of  the  parameters  for  range 
correctlor;  and  their  residual  errors  .are  listed  at  tiie  local  times  of  the  dally  moan 
maxima  and  minima  for  the  periods  representing  the  seasonal  maxima  and  minima. 
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Table  1.  Residual  Range  Error 


winter 

V ernal 
Equinox 

Summer 

Autumnal 

Equinox 

Max. 

Min. 

Max. 

Min. 

Max. 

Min. 

Max, 

Min. 

Solar  M 

aximum 

(R  = 110 
z 

r ' ' 

5’=  155) 

AR 

75. 

0 

12. 

5 

95.0 

17.5 

55.0 

17.  5 

90.  0 

15.0 

6H 

m 

13. 

5 

4. 

0 

16.  0 

5,5 

9.0 

5.  0 

16.  0 

5.0 

">3h 

9. 

0 

3. 

5 

iO.  0 

3.5 

6.0 

3.  5 

9.  0 

4.0 

'"ih 

5. 

0 

2. 

3 

5.5 

2.3 

• 1 

2.  5 

3.7 

1.8 

®R30„. 

3. 

5 

1. 

2 

3.5 

1.2 

2. 5 

1.  5 

5.  0 

1.  5 

Sola 

r Minimum  = 10;  S’  tt  71) 

z 

AR 

22. 

5 

3. 

5 

22.5 

3.5 

20. 0 

4.  5 

25.  0 

4.  5 

4. 

0 

1. 

5 

5.5 

1.5 

3.5 

1,0 

3.  5 

1.  5 

“'3h 

4. 

0 

1. 

5 

5.5 

1.5 

3.5 

1.  0 

4.  5 

1.5 

“Xlh 

2. 

5 

0. 

8 

2.5 

0.5 

1. 5 

0.  8 

2.  5 

0.  8 

'"som 

4 

2. 

0 

0. 

1.5 

0.3 

1.0 

0.  5 

1.5 

0.  5 

AR  -the  ionospheric  component  of  range  measurement  in 
feet  for  a 1300  MHz  radar  with  target  at  5“  elevation 
and  1000  km  altitude. 

5 -the  residual  error  in  range  correction  using  a median 

prediction,  caused  by  the  day-to-day  variability  of 
the  ionosphere. 

^^3h  residual  error  in  range  correction  using  a scaled 

median  prediction  3 hours  after  updating. 

-the  residual  error  in  i-ange  correction  using  a scaled 
median  prediction  1 hour  after  updating. 

RjOm  residual  error  in  range  correction  using  a scaled 

n>edian  prediction  30  minutes  after  updating. 


To  estimate  the  range  correction  for  other  values  of  S’,  it  is  possible  to  assume  a 
linear  relationship  between  S’  and  TEC  and  interpolate  or  estrapolato  for  the  desired 
value.  ^ 

These  results  show  that  the  upd.atlng  technique  based  on  the  pulse -pair  measure- 
ments car.  BuccesafuUy  reduce  the  errors  caused  by  the  day-to-day  variability  of 
the  ionosphere  .•‘or  up  :o  3 hours  after  the  initial  calibration. 


14 


•'V-- 


& COMPARISON  OF  AFGWC  AND  CONTRACTOR  PREDICTIONS 

Before  30  September  1977,  the  refraction  corrections  for  COBRA  DANE  were 
aerived  from  an  independent  model  developed  by  the  contractor  for  use  during 
initial  tests.  Using  a set  of  prediction  cards  kindly  supplied  by  the  contractor,  a 
comparison  to  the  AFGWC  predictions  was  made  for  the  month  of  September  1977. 
The  two  programs  both  depended  on  the  standard  ITS -7  8 coefficients  for  determining 
ionospheric  parameters.  The  essential  difference  in  these  two  programs  was  their 
use  of  separate  techniques  for  determining  the  height  of  the  Fg  region,  and  slight 
differences  in  the  shape  of  their  electron -density  profiles. 

A Plot  of  C^  (Figure  2)  representing  the  monthly  median  for  September  1977 
was  made  with  AFGL's  model  (sunspot  number  equal  to  20)  and  compared  to  Ray- 
thc.  n's  ...odel  (SSN=20).  A variation  of  about  5 feet  was  found  to  exist  between  the 
;.wo  lines  of  data,  with  AFGL  data  being  displaced  above  Raytheon  data,  indicating 
some  systematic  difference.  Also,  some  non-systematic  difference  was  indicated, 
as  Raytheon  data  was  higher  than  AFGL  data,  from  1900  hours  to  2000  hours  (UT). 


Figure  2.  Comparison  of  Monthly 
Median  Correction  as  Calculated 
by  Two  Models 


In  an  attempt  to  compensate  fo*-  \he  systematic  difference  fo  .kI  using  SSN=20, 
the  sunspot  number  of  AFGL  was  set  at  1 1.  Upon  compa  ison  of  SSN=  1 1 with  the 
original  Raytheon  data,  it  w.a8  found  that  this  pio'  of  GWC/AFGL  was  displaced 
below  the  Raytheon  Une  throughout  the  day.  Visual  interpretation  revealed  varia- 
tion in  excess  of  the  previous  5 -foot  variation,  especially  at  the  begimilng  of  the  day. 
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Plots  of  the  azimuth  variation  (Figure  3)  were  done  for  five  specific  time 
periods  (02,  04,  06,  14,  and  22  UT).  Data  from  the  two  models  compared  most 
favorably  during  the  daytime  (hours  02,  04,  and  06  UT).  These  points,  covering 
33  azimuth  sectors,  showed  remarkable  similarity  between  the  two  models.  The 
lines  tended  to  separate  slightly  (mostly  at  14  UT)  at  either  end  of  these  graphs. 

In  general,  Raytheon  was  displaced  above  AFGL  (especially  at  the  beginning  of  each 
line),  with  the  exception  of  22  UT  wherein  Raytheon  remained  below  the  AFGL  line 
throughout  the  plot. 

A comparison  of  the  height  vector  (Figure  4)  was  done  for  the  following  four 
time  periods:  02,  08,  14,  and  22  UT,  and  twin  plots  of  height  versus  33  values  of 
height  were  made.  In  general,  the  Raytneon  line  was  displaced  above  the  GWC  line, 
with  the  exception  of  the  first  11  readings  (of  height)  for  the  14  UT  plot.  The  least 
variation  occurred  on  the  14  UT  plot.  The  most  variation  occurred  on  the  two 
extreme  time  periods  of  02  UT  and  22  UT. 

The  prediction  for  the  maximum  correction  near  the  horizon  (Figure  3)  and 
the  variations  with  azimuth  (Figure  3)  were  very  close.  This  resulted  probably 
because  both  measurements  depended  primarily  on  the  value  of  FOFg,  which  was 
nearly  Identical  in  each  model.  The  prediction  for  the  variation  in  height  (Figure  4) 
was  quiet  different.  This  difference  was  probably  due  to  the  fact  that  height  of  the 
Fg  region  and  the  profile  shape  were  determined  Independently  in  these  two  models. 
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